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I show that massive-particle dynamics can be simulated by a weak, external perturbation on a 
potential flow in an ideal fluid. The perturbation deflning a particle is dictated in a small (spherical) 
region that is otherwise free to roam in the fluid. Here I take it as an external potential that couples 
to the fluid density or as a rigid distribution of sources with vanishing total out-flux. The effective 
Lagrangian for such particles is shown to be of the form mc^£(U'^ /c^), where U is the velocity of 
the particle relative to the fluid and c the speed of sound. This can serve as a model for emergent 
relativistic inertia a la Mach's principle with m playing the role of inertial mass, and also of analog 
gravity where m is also the passive gravitational mass. The mass m depends on the particle type and 
intrinsic structure (and on position if the background density is not constant), while £ is universal: 
For D dimensional particles i oc F{1, 1/2; D/2; U"^ /c?) {F is the hypergeometric function). These 
particles have the following interesting dynamics: Particles fall in the same way in the analog 
gravitational field mimicked by the flow, independent of their internal structure, thus satisfying the 
weak equivalence principle. For D < 5 they all have a relativistic limit with the acquired energy 
and momentum diverging as f/ ^ c. For D < 7 the null geodesies of the standard acoustic metric 
solve our equation of motion. Interestingly, for D — 4 the dynamics is very nearly Lorentzian: 
£ oc —mc'^-y~^\{'y) (up to a constant), with A = (1 -f 7"""^)"^ varying between 1/2 to 1 (7 is the 
"Lorentz factor" for the particle velocity relative to the fluid). The particles can be said to follow 
the geodesies of a generalized acoustic metric of a Finslerian type that shares the null geodesies 
with the standard acoustic metric. In vortex geometries, the ergosphere is automatically the static 
limit. As in the real world, in "black hole" geometries circular orbits do not exist below a certain 
radius that occurs outside the horizon. There is a natural deflnition of antiparticles; and I describe a 
mock particle vacuum in whose context one can discuss, e.g., particle Hawking radiation near event 
horizons. 

PACS numbers: 04.20.-q 47.10.-g 



I. INTRODUCTION 

It is well documented that the propagation of acous- 
tic waves in inviscid, barotropic, irrotational background 
flows bears some enlightening resemblances to propaga- 
tion of Ught in curved space times (see the seminal paper 
of Unruh[l] and many subsequent expansions; e.g., [2, 3], 
and the recent extensive review in [4]): The flow poten- 
tial, ?7, describing weak acoustic waves moving on a given 
background I?-dimensional flow satisfies the wave equa- 
tion 

□77^(-g)-i/2[(_g)i/VryJ.„ = 0, (1) 

where y*'" is the inverse of the matrix = 
(p/c)2/(^-i)g,„ with goo = -(c^ - v^), go. = g.o = 
— Wi, Qij — Sij, with p the background flow density, 
V its velocity, and c the local speed of sound (g = 
_[p(£'+i)/^2]2/(D-i) -g ^j^g determinant of The ma- 
trix is called the acoustic space-time metric because 
eq.(l) is identical to the massless-scalar wave equation 
in a curved space time described by the metric g^„, □ 
being the covariant d'Alembertian. This setup is used to 
simulate the propagation of light in gravitational fields. 
The analogy is, however, anything but complete. For ex- 
ample, a coordinate transformation of an acoustic metric 
does not take us to another acoustic metric. And, in the 
flat-space-time analog (homogeneous background flow at 
rest) there is no parallel with observer independence of 



the speed of light. The situation is more akin to propa- 
gation of light in the old aether. Also, there is not yet 
an analog of the Einstein equations whereby the effective 
geometry is determined by its sources. Still, the analogy, 
where it does exist, is very useful and captures some cru- 
cial aspects of photon propagation in curved space times. 
For example, it elucidates the behavior of light near event 
horizons in "black hole" geometries. It can also model 
Lorentz invariance breakdown in light propagation, etc. 
(see, e.g., [2, 5, 6]). 

Here I propose to extend this analogy to massive parti- 
cles, also in the context of ideal fluids, in the usual hope 
that it might teach us about the real processes it rep- 
resents. The interest in such models might be twofold. 
First, they provide models for emergent relativistic in- 
ertia: Starting with objects that have negligible inertia 
of their own, their interaction with the fluid puts a cost 
on their motion by endowing them with an effective ki- 
netic action and thus with pseudo-energy and pseudo- 
momentum. Such models might shed light on the origins 
of real inertia, in the spirit of Mach's principle. (I mean 
here Mach's principle in the extended sense that inertia is 
not an innate attribute of bodies but emerges as a result 
of their interaction with some omnipresent agent, such 
as a field, the vacuum, the fluid in our case, or, as in 
the original view, the totality of other bodies in the uni- 
verse.) They also permit us to study possible mechanisms 
for breakdown of the standard Lorentzian kinematics at 
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high Lorentz factors. Wc can also study possible depar- 
tures from standard inertia when the global setup of the 
fluid is changed to mimic real inertia in the context of 
the non-trivial cosmology of our universe. In fact, my 
own interest in the subject stemmed originally from the 
wish to construct mechanical models for modified inertia 
that will mimic the behavior of MOND, a theory that I 
proposed to replace the need for dark matter in galactic 
systems (e.g. [7 9]). 

Second, these models extend the usefulness of acous- 
tic analogs of light propagation in curved space times to 
that of massive particles in gravitational fields. Interest- 
ingly, I find that for properly defined particles the same 
attribute that plays the role of inertial mass also plays 
the role of passive gravitational mass thus conforming to 
the weak equivalence principle. With these models we 
can study mechanisms for the breakdown of the weak 
equivalence principle, dynamics of massive particles near 
black hole analogs, such as the existence of a last stable 
orbit, etc.. And, with an appropriate definition of the 
particle vacuum we may be able to study Hawking radi- 
ation and other quantum eff'ects in curved space time for 
massive particles. 

The particles I shall describe do not generally fol- 
low geodesies of the acoustic metric itself, for which the 
proper time is 

dr = a^-'^dt, (2) 

where 7 = {1 — [dx/dt — i/(r)]^}~^/^ is the "Lorentz 
factor" of the velocity relative to the fluid, and a = 
[pc^^~'^^]^^^^~^\ In the real world Lorentz invariance 
dictates the above path length as the particle action; 
but this is not so in the fluid context. Nevertheless, it 
would still be useful to find analogs that have enough of 
the properties of real particles, in particular, relativistic, 
quasi-Lorentzian dynamics. This I begin to do in this pa- 
per. As I shall show, there is, in fact, a generalization of 
the acoustic metric in the form of a Finslerian one whose 
path length is the particle action and which shares the 
null geodesies with the acoustic metric; so, a unified de- 
scription of massive and massless particles does emerge 
with a Finslerian acoustic metric. 

I am not concerned here with the practicability of lab- 
oratory construction of such analogs. I view their useful- 
ness mainly as theoretical laboratories for testing ideas 
concerning inertia and gravity. 

A rather different approach towards mimicking massive 
particles in the context of Bose-Einstein condensates is 
described in [10, 11]. The acquisition of induced mass by 
vortices moving in superfluids has been discussed in [12] 
and references therein. 

In section II I discuss the general idea and define the 
particles. Section III contains the derivation of the ef- 
fective particle action in D dimensions. In section IV I 
discuss various aspects of the resulting dynamics of the 
particles, first in flat space times, then in the presence 
of analog gravity. Section V brings up some additional 
issues. 



II. MASSIVE PARTICLES IN FLAT SPACE 
TIME 

Very weak perturbations of the fluid flow itself, to wit 
acoustic waves, are the analog of light in the fluid context. 
They are described by the same degrees of freedom as the 
unperturbed, background flow and move with the local 
speed of sound relative to the fluid. Analogs of massive 
particles should be able to move at any "subluminal" 
speed relative to the fluid. They should thus be deflned 
as regions of space where the equations of motion for the 
background flow break down. The exact definition of the 
particles, with the prerequisites they have to satisfy, is 
best presented in the context of "flat" space times; i.e., 
homogeneous fluids at rest. 

A rigid body moving with constant speed in an invis- 
cid, incompressible fluid is subject to no force; this is 
known as the d'Alcmbcrt paradox. It follows straight- 
forwardly from scaling that the energy of the fluid in 
this configuration is proportional to the body's velocity 
squared. Accelerating the body increases the energy of 
the fluid hence creating an effective force resisting accel- 
eration. The force is proportional to the acceleration, 
giving rise to acquired Newtonian inertia: The speed 
of sound in incompressible fluids is formally inflnite, so 
there is no velocity parameter to represent the speed of 
light. The effective inertial mass tensor is the fluid den- 
sity times some geometrical volume matrix that depends 
on the shape of the body (See e.g. [13] and [14] §11). My 
aim here is to extend this idea to bodies with relativistic 
kinematics. Naturally one then begins with a compress- 
ible background fluid whose speed of sound will play the 
role of the speed of light. A rigid body is however not 
a good model for a particle (see below) so we'll have to 
flnd others. 

A more general discussion of forces on static bodies in 
a class of nonlinear media, of which our simple flow is an 
example, can be found in [15] where I consider different 
possible deflnitions of bodies. For example, to deflne a 
body one can dictate boundary conditions on a closed 
rigid surface of a region that can move in the fluid. Dic- 
tating a vanishing normal component of the flow velocity 
defines a rigid body, for instance. Alternatively, we can 
define the particle as a rigid collection of sources (and 
sinks). Yet another way is to dictate inside the parti- 
cle an external potential that couples to the flow density. 
And yet another is to take the particle as a small region 
of non-vanishing vorticity. There are more options; the 
choice, however, is limited by the following requirements 
that I think should apply in the quest of analog massive 
particles: a. The particle should constitute a control- 
lably weak perturbation on the background flow. This 
is not just to facilitate the derivation of the particle's 
dynamics but mainly to prevent the particle from prob- 
ing the equation of state of the fluid at densities other 
than the background value. This requirement eliminates, 
for example, rigid bodies as candidate particles because 
when such bodies move with relativistic speeds they ere- 
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ate strong perturbations in their vicinity, no matter how 
small in size they arc. b. An effective principle of in- 
ertia should hold: when the particle is set in motion in 
a homogeneous fluid it should retain a constant speed; 
i.e. it should not be subject to forces by the fluid. The 
d'Alembert paradox insures this for a rigid body in an in- 
compressible ideal fluid. It was shown in [15] that this is 
true also for compressible fluids, and, in such a case, also 
when the body is defined as a region of dictated exter- 
nal potential or as a distribution of sources, provided the 
integrated source out-flux vanishes (this holds exactly, 
even when the body is not a weak perturbation), c. To 
limit the scope of the discussion I also require in this pa- 
per that the particle is a rigid object, with no internal 
degrees of freedom. It may be interesting to relax this 
assumption in various ways. 

I shall indeed concentrate here on particles defined as 
either a distribution of sources or a dictated potential. 
The former is epitomized by a source-sink dipole such 
as a small pipe within which there is a pump sucking 
fluid at one end and ejecting it at the other, or by any 
arrangement of dipoles such as a dipole source layer, etc.. 
The second type may be realized archetypically as a set of 
electric charges held together rigidly by a structure that 
does not disturb the fluid mechanically (a rigid cage), and 
moving in a weakly charged fluid with constant charge- 
to-mass density ratio. It is best to take the total body 
charge as zero so as to attain a confined potential. Ideally 
it would be good to add an inert (static) background with 
the opposite charge to cancel that of the unperturbed 
fluid, so that only density perturbations carry net charge. 



III. THE EFFECTIVE PARTICLE ACTION 

One can get the equations of motion of an irrotational, 
inviscid, barotropic fluid in D space dimensions from the 
effective action 

5 = - y d'^+^xipc/),, + ^p(V0)^ + e + pe + pij + 4>s],{3) 

where cP^^x stands for dtd^r, p is the fliud density, (j) is 
the velocity potential: v = V4>, and e{p) is the intrinsic 
energy per unit volume, which is a function of p for a 
barotropic fluid. The action S is based on that derived 
in [16] to which I have added a source term with source 
density s{f,t), and potential terms. The potential fields 
0{f, t) and 'tp{f, t) couple to the fluid density. They can be 
of the same type but I write them separately because they 
have different roles: 9 is completely dictated externally, 
and partakes in establishing the unperturbed background 
flow, while V represents a particle and so constitutes a 
weak perturbation confined to a very small, freely moving 
region of space. (I kept here the sign of the action, which 
is derived in [16] from the fluid action J pv^ /2 — e; so 
it is clear with which sign to add actions for additional 
degrees of freedom. Kinetic energies are added with a 



plus sign, while potential energies, such as p^, appear 
with a minus sign.) 

Varying the action over gives the continuity equation: 

p,t + V ■ {pVcj)) = s. (4) 

Varying over p gives the Bernoulli equation: 

^.. + l{v^f + h{p) + e + ij = o, (5) 

which in the barotropic irrotational case is equivalent to 
the Euler equation. Here h{p) = e'{p) is the specific 
enthalpy. 

To the fluid degrees of freedom we now add those of the 
model particle: its position r*(f), and possibly its orienta- 
tion. At this stage I want to eliminate the orientation as 
an unnecessary (but possibly interesting) complication. 
Later on I shall assume a spherically symmetric particle 
for which this is not an issue. For the time being I shall 
take an arbitrarily shaped particle but assume that its 
orientation is kept fixed in space (e.g. by providing it 
with a gyroscope), which generically gives anisotropic in- 
ertia. The particle's dynamics will turn out to depend on 
its orientation with respect to its velocity relative to the 
fluid. If we can somehow keep this orientation fixed (e.g. 
by providing the body with efficacious fins) dynamics will 
be isotropic for any body shape. 

One type of particle I treat is a small region of space, 
positioned around r, (i) , where a rigid arrangement of 
sources is dictated. This corresponds to a source distri- 
bution 

s{r,t) = s[r-u{t)]. (6) 

s vanishes everywhere except in a volume of diameter 
a much smaller than any length scale characterizing the 
unperturbed flow and the trajectory of the particle, and 
corresponds to a vanishing total outflow: 

j s{f)d^r = 0. (7) 

A different type of particle may be represented by a small 
volume of diameter a around r, (t) where an external po- 
tential is dictated: 

i;{f,t)=i>[r-f4t)] (8) 

with ip vanishing rapidly beyond the radius of the source. 
While s{r) or x(;{-r) are fixed and constitute the internal 
structure of the particle, its position r* is free. My aim 
is to derive an effective action for r* by solving for p 
and (j) for a given trajectory r*(f), then substitute these 
back in the action to get a functional of r* (t) that is the 
required effective action of the particle. I do this under 
the assumption that the particle is a weak perturbation 
on the background flow and is very small. 

Let us keep the designation p, (j) for the unperturbed, 
background flow attributes, and write the density and 



4 



velocity potential in the presence of such perturbations 
as p + C &iid (p + r], respectively. Expanding the action 
to second order in C and t] we get a zeroth order term, 
which is taken as a constant for a given background. The 
first order term is, after some integrations by parts, 

Si = -j d''+'x{Q[<l>, + \{Vcl>f + h{p) + 6] 
- r][p,, + V • (pV<^)] + (pr?),, + V • (r?pV(/))} 
+ j cP+^x{p'4) + (l)s). (9) 

The first two terms vanish for solutions of the unper- 
turbed field equations. The next two terms are the usual 
integrals of complete derivatives; they vanish if we can ne- 
glect the perturbation at space and time infinities. The 
second integral may engender first order effective forces 
on our particle and I want to eliminate it. For a potential 
particle this can be done by assuming a background flow 
of constant density, as I shall eventually assume anyway. 
In this case this term becomes an immaterial constant 
contribution to the Lagrangian (x p J dPrij). If we want 
to permit a variable density we add to the background 
flow an incTt background distribution of charges that can- 
cels that of p, then ijj couples only to the perturbation C, 
and this first order term disappears. For a source parti- 
cle this first order term is analogous to the energy of an 
electric charge distribution s in a potential field (j). For a 
spherical particle, such as I treat most of the time, this 
term actually vanishes for a constant background density. 
In this case the continuity equation is the Laplace equa- 
tion for (f>. Then write s cx A$, integrate twice by parts, 
and note that because s is spherical and of vanishing in- 
tegral. Gauss's theorem says that $ vanishes outside the 
source. (Our charge distribution produces no field of its 
own outside it. The external field, which satisfies the 
Laplace equation, has no sources inside the charge dis- 
tribution, so the source does not interact with the field.) 
If the particle is not spherical, expand the space inte- 
gral in this term in multipoles about r*. The monopole 
contribution 4'{r») J s{r)d^r vanishes. The dipole term 
w(r*) • / s{r)rd^r is of the same order in the particle 
size as the terms I shall want to keep (see below) and 
I get rid of it by taking particles with vanishing dipole. 
The quadrupole ,, (fv) / s{r)rirjd^r, and higher multi- 
poles, is of higher order in the particle size (while being 
of a lower order in the perturbation) and I neglect such 
contribution to the action because after the strength of 
the perturbation is set we can take the particle size as 
small as we wish. 

Turn now to the second order action, from which the 
effective particle action is constructed. We have 

Sii= - J d'^+'xicn, + ^p{Vvf 

+ CW-Vr?+|^C'+CV' + »7s], (10) 



with c the speed of sound: = p'(p) = ph'{p)). It gives 
the first order Bernoulli equation by varying over C,: 

??,* + (cVp)C + W-Vr? + V' = 0, (11) 
and the continuity equation, 

C+V-(pVr/)+V-(CV0) = s, (12) 
by varying over -q. The first can be used to eliminate 

C = -(p/c2)(,7,. + V0.Vr? + V), (13) 
and substituting in the second we get 

- [{p/c^){'n,t + W • Vry V)],. + V • {pVri) 

- V-[(p/c2)(r?,, +V.^.V?? + ^)V(/>]=s. (14) 

Rearranging gives 

- [(p/c2)(7?.*+V</)-Vr;)],,+V-(/9Vr?) 

- V-[(p/c2)(»7,.+W-V»7)W] 

= s + [(p/c2)V]., + V • [(p/c2)V>W]- (15) 

I use the continuity equation for the background flow to 
write s + p{'ijj/c^) t 4- pV0 • V{ip/c'^) for the right hand 
side. The left hand side is {—gY^'^Drj, where, as in eq.(l), 
□ is the covariant d'Alembertian corresponding to the 
acoustic metric gf^„. So 



□r? = s + (^/c2),,J^ 



(16) 



where s = {—g) ^/"^s is the covariant source density, and 
the current 

jM^(_5)-i/2^(l,^T) (17) 

is covariantly conserved: 

^^. = (-5)-'/'[(-9)'/'Jl,.=0. (18) 

Now eliminate C, from the action 5// itself. After some 
algebra, one gets inside the integral (up to a total deriva- 
tive) 



j d''+'x{-gY'^[\n,,n^^g-'' + n{'^/^),,J^ 

(19) 



Varying over rj gives the field equation (16). 

The program is then as follows: for a given r»(t), 
which, together with the given s or ij), determines the 
source term, solve eq.(16) for ri{r,t), then substitute it 
in the expression for Su to get the value of the effective 
action as a functional of the trajectory; this is the par- 
ticle action we are after, S'[7v(i)]. Equation. (19) requires 
knowledge of the solution rj everywhere in space time. A 
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more manageable expression is gotten by employing the 
integral relation 

= 0, (20) 

which holds for solutions of the field equation up to sur- 
face terms at space-time infinity. (This is a simple special 
case of the results of [17] and follows straightforwardly 
by integrating the first term by parts to give —{l/2)rj\3r] 
then using the field equation (16).) So we can set 



2pc^ ^ 



(21) 



This expression requires knowledge of rj only inside the 
particle, where either ijj or s don't vanish; this is very 
helpful. 

It is impracticable to solve for rj for an arbitrary tra- 
jectory in an arbitrary background fiow. It is clear that 
the resulting eff'ective action would be time non local. 
However, the assumed smallncss of the particle permits 
us to approximate ij inside the particle in a way that de- 
pends only on the instantaneous state of motion, and this 
will result in a local approximation of the action. For a 
very small particle we can assume that as it moves about, 
a steady state corresponding to the instantaneous condi- 
tions is reestablished within the particle on the short time 
scale it takes sound waves to get from one end of it to the 
other. We essentially separate the dependence on macro- 
scopic coordinates and the microscopic ones within the 
body, where r] changes quickly, by assuming that from 
eqs.(6) and (8) we can write to a very good approxima- 
tion (becoming exact in the limit of infinitesimal particle 
size) 



r]{r,t) = fi[r- r^{t)] 



(22) 



to describe the fast variations of rj around the parti- 
cle's position in space time, and where fj still depends 
on macroscopic properties such as the flow and particle 
velocities and the fluid density at r*(t). I shall discuss 
below the conditions for this approximation to hold. 

We now calculate Or] with this ansatz. I again make 
use of the fact that, due to the smallness of the particle, 
the space and time variations of rj are dominated by those 
produced by the fast variations in f/[r — r*{t)]). So, for 
example, in rj t = —v^, ■ Wf] + q, where q represents terms 
coming from the implicit dependence of fj on macroscopic 
quantities and their time variation, we neglect all such 
terms. Then in rj^t.t ~ {v* • V)^f} — (dv^/dt) ■ Vr) (again 
neglecting q terms) I further neglect the second term (by 
our approximation \dvt:/dt\ <C v^/a generically, a being 
the diameter of the particle). With this approximation, 
which leaves us only with terms with second derivatives 
of 77, we get 



□r?=^{A-[^.Vm 
P c 



(23) 



where, U = v*—vis the relative velocity of the particle 
with respect to the fluid, and v, p, c are evaluated at 
n{t). Thus 

fj,.,. + fj,,,, + 7"\„. = p-^s - (f//c2)V;,., (24) 

where 7 is the relative "Lorentz factor" 7 = (1 — 

?7^/c^)^^/2, the z axis is in the direction of U, 
and where I used the fact that now J'^{tp/c'^)^^ = 

(-fi')"^^^(p/c^)(V',t + vVip) = -{pc)-^U ■ yi>. In the 
coordinates x' = x, y' = y, z' = 'yz eq.(24) becomes the 
Poisson equation 



A't) = p-^s[R{R')] - Uc-'^4>JR{R% 



(25) 



provided the relative speed U is subsonic. This appear- 
ance of the stretched Laplacian in the linearized equation 
for weak perturbations moving with subsonic speed in a 
compressible fluid is familiar, for example, from the treat- 
ment of a constant flow past a thin wing very nearly par- 
allel to the flow (e.g. [14] §124). The perturbation there 
enters not through source terms as here, but through the 
boundary conditions on the rigid wing, leaving us with a 
distorted Laplace equation instead of Poisson's as here. 
This equation is elliptical for subsonic speeds for which 
our treatment below applies, but become hyperbolic for 
supersonic speeds. 

The effective particle action can then be written as 



S[r*it)] = Jldt, (26) 
with the particle Lagrangian 

L = - ^Jd^'r 7?(r)s(r-) + J d^'r f}{U ■ V)^, 

d^r (27) 



+ 



2c2 



where (for D > 2) 



-2)nD J 



s — c 



R 



•i\D-2 



(28) 



is the solution of the stretched Poisson equation; is 
the solid angle in D dimensions. The D = 1,2 cases will 
be treated separately in Appendix A. 

To recapitulate, the approximation I made amounts to 
the following procedure: At any given time take the local 
values of the velocities of the fluid and the particle and 
of the fluid density, calculate the steady state solution, 77, 
from eq.(24) for a homogeneous fluid with these proper- 
ties and an eternally constant particle velocity, then use 
this for the instantaneous r] inside the particle. 

I now proceed to discuss separately source and poten- 
tial particles. 
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A. Source particles 

For a pure source particle put = 0; then substituting 
expression (28) in eq.(27) and changing to the f* variables 
we get (for D > 2) 



L = 



1 



2{D - 2)QDn 



J d'^r'd 



-• (29) 



can be obtained analytically. This is done in Appendix 
A and yields 



L = ioF(l, 



1 D C/2 
2' 2"'"^ 



(34) 



where Lq = Eo/i^oP is the value of the effective La- 
grangian for J7 = 0, and F is the Gauss hypergeometric 
function. 



The integral in eq.(29) is proportional to the "electro- 
static" energy of a charge distribution s stretched by a 
factor 7 in the z direction. Note that L is positive be- 
cause it is proportional to — / ryAf) = J (Vf))^ (^'^ is a 
positive definite operator). 

We can also write the integral in terms of the f coor- 
dinates as 
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2(£>-2)OdP 

d^rd^R s{R)s{r} 
{{r- RY + 72[([//c) • (r - i?)]2}(i?-2)/^ 



,(30) 



where the full dependence of L on 7, the structure of the 
particle, and its orientation with respect to the relative 
velocity is explicit. 

In the non-relativistic limit, C/ -C c, eq.(30) tells us 
that 



L = 



-^ + ic/,m,,C/,+0(C/V); (31) 



Eo = [2(D - 2)]-^ Jd^rd^R s{R)s{r)/\r- 
is the value for the unstretched configuration, and the 
effective mass tensor is 



B. Potential particles 

Consider now a pure potential particle (s = 0). Re- 
peating the same argumentation as before 



L = 



2{D - 2)f2Dc47 



1? -i?'P-2 



+ 



(35) 



Or, with derivatives with respect to f and R', 



L 



PI 



2{D - 2)Or,c4 



+ 



I 



IP -i?'P-2 



(36) 



In the non-relativistic limit 



ma = 



2{D - 2)nDpc^ 



I 



|f-i?P-2 



(r — Rj-' 

In the isotropic case (for which a cubic symmetry of 
the particle suffices) we get the mass of the particle 



m = Tr{mij)/D = 



2Eo 



Dflopc 



,2 ■ 



(33) 



When the background density is not a constant of the 

configuration this mass parameter is a function of space- 
time position through p and possibly c. I shall still refer 
to it as the mass of the particle. 

To insure isotropy of inertia I shall assume henceforth 
that our particle is spherically symmetric. In this case L 



L=\^j ^^d'^r + ^U,m,,U, + 0{U^/c% (37) 
where the mass tensor is (integrating by parts) 



/ 



{D - 2)Q.DC^ 



d^r d'^RiP{R)i}{f)^^\f-R\-^^-^\ {38) 



dRi drj 



In the isotropic case 



m : 



(39) 



where I used -A[r-(^-2)] = {D - 2)f]£,(5^(f). The first 
integral in eq.(36) may be viewed as proportional to the 
energy of a polarized medium with unidirectional polar- 
ization of magnitude oc U'tp[r{-r')]. It can be calculated 
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exactly for an arbitrary, spherically symmetric distribu- 
tion ipii")- The integral is calculated in Appendix B, and 
when added to the second integral we get 

= ^mc'Fill^;^), (40) 

with the same hypergeometric function appearing in the 
Lagraiigian of a source particle. These identical results 
arc obtained from rather different starting expressions, 
and I have not been able to find an underlying physical 
reason for the equality. 

C. Aspherical and compound pEirticles 

A larger variety of 7 dependences of L is afforded by 
considering aspherical particles. If the orientation of the 
particle is kept fixed in space, anisotropic inertia results 
generally; but, if we can somehow keep the orientation 
fixed with respect to U, the effective inertia is isotropic. 

As an example, consider a hyper-planar, bipolar source 
layer (charged planar capacitor in the electrostatic ana- 
log) whose normal always makes an angle Q with the 
relative velocity vector. (Actually, because of the reserva- 
tions discussed above, we need to take two, back-to-back 
dipolc layers to annihilate the dipole moment of the parti- 
cle, but this is immaterial for the results since the two lay- 
ers do not interact, so I shall just continue to speak of one 
layer.) From eq.(29) the Lagrangian is 'y~^Ec, Ec being 
the energy of the stretched bilayer. This energy, like that 
of a charged capacitor, is Ec oc Q'^d/A, where Q is the 
total charge on one layer, A the area, and d the spacing. 
Under stretching Q 7(5, A A(cos^6 -|-7^sm^9)^/^, 
and d dj{cos'^Q + 'y^sin'^Q)-'^/^. So, 

Vos2e + 72sm2e i-{u^/d^)cos^e' ^ ' 

This varies between L = Lqj^ for 6 = 0, as in the 1-D 
case, to a constant L = Lq when the bilayer moves par- 
allel to itself relative to the fluid. Integrating over angles 
with weight sin^~'^QdQ gives back our result for the 
spherical case. (In the relativistic limit the contribution 
to the action of a spherical particle then comes from a 
small © region near the leading point on the sphere. The 
corresponding area decreases with increasing dimension; 
hence the strong D dependence of the relativistic limit.) 
All the above applies to any collections of bi-layers mak- 
ing the same angle with the relative velocity vector, for 
example a cone of half-opening angle 7r/2 — G moving al- 
ways along its axis relative to the fluid, like an arrowhead. 
In general, if we tie together several particles that are so 
far from each other that their mutual interactions can be 
neglected compared to their self interactions, the body 
will have an effective Lagrangian that is the sum of those 



of the individual components. (The inter-component dis- 
tance still has to be small compared with the scale over 
which macroscopic properties of the flow vary). 

D. Limitations and caveats 

I made two types of assumptions about the parti- 
cle: that it constitute a very weak perturbation on the 
fluid even within the particle itself, and that it is very 
small in size so that conditions in it very quickly take 
up steady state values corresponding to the momentary 
macroscopic conditions around it. I now discuss in more 
detail what these require from the strength of the poten- 
tial tp, from the source density ,s, and from the particle 
size a. 1 find that in general these approximations break 
down at high 7 however small a and •ip or s are. The per- 
turbation treatment assumes that ( <^ p. Equation(13) 
gives 

C = -{p/c'Kv.t+v-yrj + ij) 

« (p/c2)(c7.V7?-V'), (42) 

where I used our approximation in the second equality. 
Consider first a potential particle. For [/ — > 0, 7? — > 
and the basic requirement is \C\/p « ^ 1- We 

further have to insure that for the limit U c wc still 
have (/p ^ 1, so we need \U ■ Vryj/c^ ^ 1 everywhere 
in the body. The exact constraint this puts on 7 de- 
pends on the particle structure. We can get an estimate 
of this quantity by noting that what I calculated as the 
first term contributing to L is {p/2(?) J d^rfj{U ■ V)'4> = 
— {p/2c^) J d^r'ijj{lJ-V)fj. If we take, for example, a parti- 
cle of constant ip we know the value of {U-V)fi is constant 
inside the particle and from the results for L it is 

[/•Vr?=-Vi[F(l,i;|;^)-l]. (43) 
So we also need 

(H/c2)[F(l,i;|; J)-l]«l. (44) 

In the non-relativistic limit the expression in s([uare 
parentheses behaves as U'^/c^, so no new requirement 
is added. In the relativistic regime we can get a validity 
limit on 7. For D > 3 the F above is finite for [/ = c 
and we do not get an additional constraint; the basic one 
suffices for all values of 7. For D = 1 we have to have 
< 1, for D = 2: j\^^\/c^ < 1, and for D = 3: 
ln{'y)\tp\/c'^ <C 1. For a source particle we have, as be- 
fore, to first order in the source strength ( w {p/c^)U-'Vfj, 
which vanishes when the particle is at rest with respect 
to the fiuid. Consider such a particle in a static fiuid. Ne- 
glecting the variation of the external potential across the 
particle, the Bernoulli equation tells us that C is indeed 
second order in the source. If we write s = pA<l> we have 
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from the continuity equation that to first order ?) = $ 
(with 3rd order corrections), and so ^ ~ — (p/2c^)(V$)^. 
The basic requirement from s for our approximation to 
hold is then |V$| ^ c. And here too, for 7 1 we 
have to have \U ■ VtjI/c^ <S 1, which casts a constraint 
on 7 that may depend on the particle structure and the 
dimension. 

Consider now the condition on the particle size. When 
the particle, having diameter a, is moving with velocity 

U relative to the fluid, it takes a sound wave time St « 
(a/c)(l-[//c)-i = (a/c)(l + [7/0)72 to move from the aft 
of the particle to its fore. We want 6t to be much shorter 
than any time scale, T, over which the environmental 
parameters change. The basic requirement, which should 
hold even at low relative velocities is then a <C cT. In the 
relativistic regime the requirement is that a ^ cT/2'y'^. 

We thus see that even for very small values of tp or s, 
and of a, our approximations, and thus our results, are 
not valid for 7 — > 00. This is to be expected: For D < 5 
the particle energy diverges for 7 ^ cxo, which, if valid 
indefinitely, says that we cannot accelerate our particles 
to supersonic speeds; but this is clearly not true. 

It would be interesting to see how our acquired dy- 
namical properties of particles arc modified when the ap- 
proximations break down, as some of these may also be 
taking place in reality. For example, the time locality 
of the Lagrangian is only a result of the approximation: 
The effects of the particle on the fluid at one time affect, 
at some level, the motion of the particle at another time. 
Indeed they may affect other particles as well, thus creat- 
ing an effective interaction between particles mediated by 
the fluid akin to the Cooper pairing interaction between 
electrons in a superconductor. We also have here some 
ready made mechanisms for the breakdown of Lorentzian 
dynamics at high 7. All these departures are still con- 
sidered in the context of inviscid, irrotational, barotropic 
fluids; and these attributes are also only approximations 
(see a discussion of these in the context of phonon prop- 
agation in [6]). 

I have also neglected the goings on inside the source 
itself. This is after all some parallel flow (e.g. in some 
pipes with pumps) that move the fluid from the sinks to 
the sources. 



IV. PARTICLE DYNAMICS 

Notwithstanding the absence of true inertia in our par- 
ticles, they acquire relativistic inertia through their in- 
teraction with the fluid: this is encapsuled in the kinetic 
action for free particles. 

From now on I shall assume a position and time inde- 
pendent density for the background flow. This situation 
is rather less cumbersome to describe, captures most of 
the concepts I want to introduce, and insures a constant 
speed of sound, which after all is our analog of the speed 
of light (this can also be insured, without imposing a con- 



stant p, by having a fluid equation of state of the form 
p = c^p + const.). It also means that m is a constant 
and so both types of particles have the same motion. (If 
the fluid density depends on position or on time we get 
a variable mass for the particles, which might be inter- 
esting to explore.) The freedom left in selecting a back- 
ground flow-a solution of the fleld equations with s = 
and tp = 0-is then only in choosing the velocity poten- 
tial fleld from among the harmonic flmctions. We then 
have to impose an external potential 9 that will satisfy 
the Bernoulli equation for the chosen velocity field. The 
exact form of the equation of state is immaterial since we 
shall probe it only at one density value where its deriva- 
tive only (speed of sound) has to be known. 

For either definition of a particle the effective action 
can be written in the form 



S = 



H[j{t)]dt. 



(45) 



For a source particle m is proportional to the "electro- 
static" energy of a charged medium with charge density 
s, and for a potential particle it is proportional to the 
"magnctostatic" energy of a sphere with unidirectional 
polarization ijj. 

I now proceed to discuss various aspects of the result- 
ing dynamics. 



A. Flat-space-time dynamics: emergent inertia 

In a flat space-time; i.e., in a homogeneous background 
flow at rest, the effective Lagrangian for the two types of 
particles discussed here is 



L = mc\D/2)F{l,^;^;^p, 



(46) 



where is the particle velocity and m is constant. Note 

that D here is determined by the symmetry of the parti- 
cle, and is not necessarily the dimension of the space in 
which it moves. For example, a plane symmetric particle 
that moves only along its normal has D = 1. More gener- 
ally, a particle in N dimensions of cylindrical symmetry 
having the symmetry of Cg) R^^^-^^ whose velocity is 
in the subspace corresponds to dimension D. 

Using the formula for the derivative of the hypergeo- 
metric function we get for the momentum 



P ■ 



dL 
9v* 



3 D + 2 



(47) 



The kinetic energy 



D, 



-L + Lo + v, ■ p = mc'i-il - F{1, ^ 



3 D + 2 vl 



2' 2' c^" 
(48) 



where I added a constant so as to make Ef. vanish for 
= 0. 



9 



These energy and momentum arc conserved as follows 
from Nother theorem's related to the assumed time inde- 
pendence and homogeneity of the fluid: If the particle is 
subject to a conservative force derived from a potential 
^, we have to add / —^[f^,{t)]dt to the action we started 
with; so; the particle now satisfies ^ = — V^, and Ek+£, 
is conserved. And if we have some inter-particle forces 
'^Pi is conserved, and if these forces are derived from a 
potential again the total energy is conserved. 

Of course, E and p are not the real energy and mo- 
mentum of the particles; these were assumed to have no 
inertia of their own so they can carry no energy and mo- 
mentum. It is the stirring of the fluid by the motion of 
the particles that puts a real inertial cost to their motion. 
The rates of change of Ek and p equal the rates of c;liangc; 
of the energy and momentum of the fluid induced when 
the particle changes its velocity; they are thus equal to 
the external power input and the external force imparted 
to the particle. Such quantities are called pseudo-energy 
and pseudo-momentum. These are often useful in the 
description of motion of objects in homogeneous media 
with which they interact (See [18] 2.4 and the review by 
[19]). A direct ciakudation of the energy and momentum 
of the fluid might have also served, but it is impracti- 
cal. Attempting to calculate them even for the simple, 
steady-state configuration, with the particle ever at con- 
stant velocity, gave me ambiguous results: When these 
quantities are written as integrals of fluid attributes over 
a volume that has to be taken to infinity, the results 
depend on the shape of the integration volume. This is 
similar to what Peierls [18] finds when trying to calculate 
these quantities for a sound wave. The present approach 
of proceeding through the action seems to be the proper 
way to proceed. 

Since F{l, l; f ;^§) = 1 + D-^v^Jc^ + 0{vt/c^), we 
have in all dimensions the non-relativistic behavior 

L-Lqk ^mvl, p w miT*, Ek w ^mvl. (49) 



For £) = 2 



L = mc 7, p = my v^, 

2 7^(7' +7-1) 2 3 

Ek =mvi > mc 7 . 



7+1 



For £) = 3 



2 3 fl + v^^/c 
L = mc — — r^ln ' 



4(t;*/c) vl - v^^/cj ' 



(51) 



(52) 



p = mv* 



1 , fl + V^:/c 

In ' 



2(t;*/c)2^' 2{v*/c) \l-v*/^ 



and 



3 . 2 

-mu*7 , 



3 2ri 1 7 A + V*/C 21 

Ek = -TOC^ 1 - — — ^ +7 

2 {v*/c) l-v*/c 

^ 2 2 

-mc 7 . 



(53) 



(54) 



The case D = 4 is particularly interesting. We can 
then write, using formula 9.131.2 in [20] (henceforth GR) 

L =2mc'F{l,^;2;^) 

= Amc'[F{l, i; 1; 7"^) - 7-'i^(l, ^; ^; 7-')].(55) 



2' 2' 



However, we have generally F{l,b;b;z) = {1 — z) ^, 
which gives 



Amc^ 
1+^' 



p = 



Amjv* 



(1+7-1)2' 



Ek = 4mc 7 



,2 (.7-2 , 1 



7+1 2 



(56) 



In the highly relativistic regime the behavior depends 
strongly on the dimension. Dimension I? = 3 is critical in 
some sense: Because F(a, 6; c; z) is finite for z = 1 when 
c > a -h 6, L is finite as 7 — > cx) for D > 3, diverges log- 
arithmically for D = 3, and diverges as a power of 7 for 

< 3. Dimension = 5 is another critical dimension 
above which the energy and momentum remain finite for 
7 ^ 00; for D = 5 itself these quantities behave as ln{y) 
in this limit (see below), while for D <b they diverge as 
a power of 7. 

Following are the Lagrangian, the momentum, and the 
energy for dimensions D < 5 in closed forms with their 
relativistic limits: 

For £) = 1 



L = mc 7 /2, p 
Ek =mvWin^ -1/2) 



TO7^W* , 



The kinematics is quasi-Lorentzian and becomes 
Lorentzian in the limit of high 7, with L « Mc^{l — y~^), 
p « Mjv^,, and Ek ~ Mc^y, where M = 4m. 

For £» = 5, using formula 9.137.14 and then 9.121.1 in 
GR we can write 

L =mc2(5/2)F(l,i;U) 

= ^^e2(.,/c)-2[l-7-2F(l,l;||)], (57) 



where 



r./i 1 3 t^*, 1 , fl+V^/c\ , 



2 4 

mc 7 . 



(50) 



2'2'c2' 2{v^/c) Vl-^^*/c 



is, in fact, the D = 3 Lagrangian. So, for 7 — > 00, p 
{15/2)mvjn{y), Ek « il5/2)mcHn{'j). 
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B. Antipeirticles and the peirticle vacuum 

If we wish to push the analogy beyond the dynamics 
of isolated, ever-existing particles, and discuss pair cre- 
ation and annihilation we need to define new notions. 
We have to identify antiparticles, and we need to have a 
proper definition of the rest mass of our particles. Also, 
unlike phonons, our particles are not an organic part of 
the fluid; they cannot be created out of it if not put in by 
hand. So, pairs will not spring OTit of the fluid even un- 
der energetically favorable conditions, such as near event 
horizons or in strong "electric" fields unless we prepare 
a "vacuum" that has the potentiality to beget them. 

For a particle given by some f/', or some s, it is natu- 
ral to define the antiparticle as that given by —ip or ~s 
respectively. They have the same m as the particles. If 
we take a particle and its antiparticle and superimpose 
them on each other we get an object that does not in- 
teract with the fluid, and in whose presence the energy 
is that of a fluid devoid of particles. The particles them- 
selves arc still there as they arc some imposed external 
structures that do not physically annihilate each other, 
they only cancel each others influence on the fluid when 
they coincide. When one separates a superimposed pair 
it increases the energy of the fluid by an amount that 
we should identify as 2moc^, with toq the rest mass of a 
particle. Let the two coalesce with no outside help, and 
the fluid goes back to its initial state plus waves carry- 
ing the released energy 2moC^ emitted as "annihilation 
radiation" . 

I thus envisage the particle vacuum as a fluid fllled with 
superimposed (annihilated) pairs. In this state they do 
not affect the fluid, and their introduction into the fluid 
does not cost in energy. The partners of each pair stick 
together, up to fluctuations, because it is energetically 
favorable. And now, if the proper circumstances arise 
pairs can be created from this vacuum and annihilated 
into it emitting phonons. 

What is then the rest mass of the particles? I have not 
been concerned so far with the energy it costs to insert 
the particles into the fluid but only with changes induced 
by their movement; so, the rest mass cannot be read off 
the Lagrangian. To determine the rest mass I considered 
in more detail a source particle whose rest mass we should 
take to be Er/c? , where Er is the energy required to in- 
troduce the source into the fluid at rest. Compare the 
energy of the steady-state fluid conflguration with the 
source inserted with that of the conflguration with an 
homogeneous fluid at rest having the same total mass. I 
argue in Appendix C that for a spherical source of van- 
ishing total out-flux Er = Eo/pQn and hence the rest 
mass is given by 



D 

mn = — m. 
2 



(59) 



Again, no new mass parameter is introduced and m also 
determines the rest mass. 



We can thus write the complete expression for the en- 
ergy of a source particle 



3 D + 2 t)2 



(60) 



C. Curved-space-time dynamics 

When the acoustic space-time is not flat i.e., when the 
background flow velocity is not constant-all our particles 
fall in the same way in a gravitational field thus obeying 
the weak equivalence principle. This is non-trivial and 
might have well been otherwise. For example, the first 
order terms that we arranged to be absent could destroy 
universal free fall. But barring such departures, which we 
saw can be avoided by properly defining the setup, m can 
be thought of as both the inertial mass and the passive 
gravitational mass of the particle. (When the background 
density is not constant the two types of particles see two 
different, but conformally related, space times.) 

Our Lagrangian is of the form L = L(U'^/c^), which 
is also true of the standard acoustic line element. These 
give the Euler-Lagrange equation 



d{L'Ui) 
Jt 



+ L'UkVk,i = Q. 



Using the fact that v is irrotational we have 
U dt dt 2 



(61) 



(62) 



Using the formula for the derivative of the hypergeomet- 
ric function we have for our actions 



L' =-mc^F{2,-; 



2' 2 



(63) 



and 



L" 



^mc^F(3,-;^;^). 



(64) 



In the limit U <^ c the first term in eq.(62) is of a 
higher order mU/c {L' and L" are finite there) and we 
are left with the standard non-relativistic equation for all 
L(?7Vc2): 



dv* ^ 
-dt = 



(65) 



where x = —v'^/2 can be identified as the Newtonian 
gravitational potential. This holds when the gravita- 
tional field is weak (v <i; c) and the motions are slow 
(f* ^ c), but also when only U <^ c. This can also 
be gotten directly from the action, which for [/ <C c is 
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/ U'^dt up to a constant. This limiting behavior is com- 
mon to all theories with an ^(7) Lagrangian including 
the standard acoustic line element. 

Because the background flow is assumed to have a con- 
stant density we see from the Bernoulli equation that 
X equals the external potential used to establish the 
background flow. (This might point the way to intro- 
duce dynamics for the acoustic metric through that of 
the external potential 6.) 

Consider now the null geodesies of the acoustic metric 
characterized hy U = c, or dr = 0. In theories for which 
L" /L' diverges when U'^/c^ 1 the equation of motion 
implies that d{U'^)/dt = if initially U ~ c, so U remains 
constant at this value. In other words, in such theories 
the null world lines of the acoustic metric are solutions of 
the equation of motion (this can be shown to hold even 
in flows with variable background density). This is the 
case for the acoustic Lagrangian itself, but also for our 
Lagrangians when D < 7. For D > 7, L, L', L" are 
finite as 7 — > 00; for 5 < D < 7, L, L' are still finite, but 
L" diverges like 7<^-"> (for D = 7 L" diverges logarith- 
mically) and so does L" /L'; for D < 5, L' diverges like 

so L"/L' behaves as 7^. 

Note that f/ = is a solution (i.e., the body just mov- 
ing with the fluid). This is also a geodesic of the acoustic 
metric. 

We saw above several instances of solutions of our field 
equations that are geodesies of the acoustic metric, and 
we shall see another in the next subsection. But in gen- 
eral the solutions of the Euler-Lagrange equations are not 
geodesies of the that metric, as the particle action is not 
its arc length. It may be useful, however, to generalize 
the acoustic proper-time interval and use our action to 
define a Finslerian one (defined only for intervals that are 
time- and null-like with respect to the acoustic metric) 
whose geodesies are the solutions of our Euler-Lagrange 
equations. For D < 3, where the Lagranian diverges for 
7 — > 00 this is well defined only for time-like elements, 
but for D > 3, df (X [£{U^/c^) - i{l)]dt is well defined 
also for null intervals of the acoustic metric (for which 
it vanishes). Furthermore we saw that for D < 7 this 
scheme embraces both massive and massless particles. 

I demonstrate this for the more interesting case = 4: 
Subtract from L a constant equal to its relativistic limit: 



(logarithmically for D = 5) and L" still as 7''' °' , 



L = L- Amc^ = -Mc'^-f-'^X{-f), 
where M = 4m, and 

A(7)=7/(l + 7). 



(66) 



(67) 



( When the background fluid density is not taken as a 
constant the term we subtract in eq.(66) is position de- 
pendent, we can then consider it as an external potential 
for the particle, which does not couple to the fluid.) We 
can then deflne a Finslerian line element 



defined for time- and light-like intervals, where here 7 
stands for {1 - [dx/dt - u(r)]2}-V2, Clearly, df{dx,dt) 
is homogeneous of order one as required. Our particles 
follow geodesies of this Finslerian metric since the action 
is S = —M(? J dr, and furthermore, phonons follow its 
null geodesies, since St = ^ dr = 0, where dr oc 
'y~^dt is the standard acoustic line element, and we saw 
above that these too extremize the Finslerian arc length 
(because they solve the equation of motion). 



D. Circular orbits in spherically symmetric space 
times 

Consider now circular orbit in a spherically symmetric, 
or axi-symmetric, configuration. Since is constant we 
can write from eq.(62) 



dv* 1^ 2 

dt 2 



(69) 



identical to the non-relativistic equation of motion. 
These orbits are thus geodesies for all choices of L{U'^/c^) 
including the acoustic one and all of our Lagrangians. For 
a circular orbit ^ = —v1r~'^r, and v"^ is a function of r 
{v is not necessarily radial). So the relation between the 
velocity and the radius is: 



(70) 



where v = \v\, f 



din 



din r ' 



Because I assumed a con- 



stant density the continuity equation and zero vorticity 

condition dictates that / is determined by the symmetry 
of the fiow. For example, in purely radial fiow in D space 
dimensions f = D — 1. 

When V is radial \U\ = v{l + /)^/^, so subsonicity dic- 
tates for massive particles that v < c(l + f)^^^^, which 
sets the limit for the innermost circular orbit. Equality 
corresponds to a "photon" orbit. In the canonical acous- 
tic black hole configuration (e.g. [6]), where v = c{rfi/r)f 
with rh the horizon radius, this implies r > rii{l + f)^^^-^ 
for massive particles. This is analogous to to the ra- 
dius occurring at 3m = 1.5rh for real Schwarzschild black 
holes in 3 -I- 1 dimensions. 

In a vortex geometry in (2+1) dimensions 



v/c-- 



-rhf/r + {rl 



(71) 



(e.g. [6]) where is the radius of the ergosphere, rh is 
that of the event horizon and eg is a unit vector in the 
azimuthal direction (I take an ingoing flow but the results 
below are the same for an outgoing one); so / = 1. The 
minimum radius of a circular orbit is 



(72) 



dT{dx, dt) = 7 ^X{'y)dt, 



(68) 



where Vj^ is for prograde motion and for retrograde 
one. 
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As in the real world is a static limit since = is 
not permitted below it lest U become supersonic. Also for 
r <rh v^, must have a component in the radial direction 
of the flow. 

I haven't fully checked the question of stability of the 
circular orbits. But note that in our constant density 
configurations the Newtonian potential X is a power law 
of the radius: x oc in the D = 2 case and with 
a higher power in higher dimensions. This means that 
the effective radial potential (including the centrifugal 
barrier) for a particle with given angular momentum has 
only a maximum for D > 2. So there aren't any stable 
bound orbits in the nonrelativistic case and the circular 
orbits are unstable. For D = 2, depending on the value 
of the angular momentum an orbit is either unbound or 
goes through the origin. This is not necessarily so in 
configurations with non-constant background densities, 
but their discussion is beyond what I wish to consider 
here. 



V. DISCUSSION 

Evidently, flow models can provide quasi-realistic 
analogs for relativistic inertia of massive particles and 
their behavior in gravitational fields. I am presenting 
these models in the twofold hope that they can inspire 
us in understanding the origins and the validity limits 
of genuine inertia, and that by considering how these 
models respond to tweaking we can learn about possi- 
ble modification of standard physics in the real world. 
For example, looking where and how our approximations 
break down we can gain insight as to where and how 
standard dynamics may go awry. Such departures may 
include breakdown of standard dynamics at high 7, time 
non-locality of the particle action, and fluid-mediated in- 
teractions between particles, all of which are not part of 
standard dynamics. Such models can also be used to en- 
lighten us on how local dynamics might be affected by 
cosmology at large. Cosmological expansion may be in- 
cluded in the context of fluid analogs and could model, 
for instance, cosmological variations of particle masses 
(through variations in the fluid density, which enters the 
normalization of the masses), or variations in the speed 
of light. My hope in this connection is to simulate the 
dynamics implied by MOND, which revolves around an 
acceleration constant, oq that turns out to be of the order 
of the cosmic acceleration. 

The fact that the more realistic models emerge for 
higher space dimensions than we seem to be living in is 
not disconcerting. Recent work on membrane universes 
has taught us that while most of the physical objects we 
deal with may be confined to sub-manifold of lower di- 
mensions some aspects of physics, such as gravity, may 
be probing the higher dimensional aspects of space-time. 
In our models we could, for example, envisage particles 
moving in a fluid in D dimensions space, but somehow 
confined to reside in a three-dimensional sub-manifold. 



This would give rise to D-dimension inertia in a lower 
dimensional effective space. 

My main purpose in this paper is to demonstrate the 
concept: Instead of considering weak perturbations that 
are part of the background itself, and which thus move 
with a speed dictated by the characteristic speeds of the 
background, define the particle as an externally dictated 
perturbation that breaks the field equation of the back- 
ground, but that can otherwise move freely on the back- 
ground field. Various extensions and generalizations sug- 
gest themselves that are worth exploring. For example, 
we can define other types of particles, or permit non-rigid 
particles with responsive intrinsic structure. This would 
produce longer range interactions between particles simi- 
lar in nature to van der Vals interactions between neutral 
charge distributions (our rigid particles interact only on 
contact). Such particles with dynamical, internal degrees 
of freedom; e.g., with the different charges in a source sys- 
tem connected by "springs" , may also serve as detectors 
for (phononic) Unruh radiation, as the internal degrees 
of freedom will couple to the phononic field. 

And, we can generalize this idea to the whole gamut of 
analog models for which photon propagation can be sim- 
ulated (e.g., [4]). One possibility, for example, is to look 
at small charge distributions (of vanishing total charge) 
in the context of non-linear electrodynamics where the 
action of the electromagnetic field is not the invariant 
Fn^F^" but some function of it. 

Finally note that as things now stand, our particles 
cannot serve as sources for a mock gravitational field 
through their effect on the fluid: Their mass m is not an 
active gravitational mass. So, they do not help towards 
constructing an analog of the Einstein-Hilbert action. 

APPENDIX A: CALCULATION OF THE 
LAGRANGIAN FOR A SOURCE PARTICLE 

Here I calculate the energy integral in cq.(29) for a 
spherical distribution of sources s(r). Divide the distri- 
bution into concentric thin shells of radii rj and total 

charges q^. The integral is twice the electrostatic energy 
of the system made of these shells all stretched by a fac- 
tor 7 in the z direction. The stretching is only of the 
geometry without thinning the density. Each spherical 
shell becomes a homocoid: a shell bound by two concen- 
tric, oriented ellipsoids of the same axes ratio 7, with the 
original density inside: so, the resulting homoeoid i has 
a total charge Qi = jqi. Write now the energy as 

i i<j 

where Ei is the self energy of shell i and Eij is the interac- 
tion energy betweens shells i and j [i is interior to j). A 
homoeoid produces a constant potential inside its cavity 
[ipi for homoeoid i) and if the homoeoid is thin, as here, 

this is also the potential on the shell. Thus, Ei = QiLpi/2 
and Eij = Qi<pj. We can most easily calculate ipi as the 
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value of the potential at the center of the cavity and this 
is simply (for D > 2) 



"D Jo 



[l + (72 - l)cOs26»](^-2)/2 



^.(1)%^27^-^ X 

/2 sin^-'^e de 



[1- ([//c)2szn26'](^-2)/2' 



(A2) 



where iln is the £)-dimensional solid angle {Q.d = 
CId-1 sin'-'~'^6d6.), and (pt{l) = qi/rf~'^ is the po- 
tential for the unstretchcd shell. The integral can be 
expressed using a Gauss hypergeometric function (using 
formula 3.681.1 in [20], henceforth GR): 



/ 

Jo 



■^'^ _ 1 n-1 l.p,D-2 D-1 D 



It can be shown that B[{D - l)/2, 1/2] = nD/O.D-i. 
Also use formula 9.131.1 in GR to further simplify and 
get 



1 D [/2, 
2' 2"'^ 

so for the Lagrangian one has 



E = Eo^F{l,-;-;^ 



L = LqF{1, -■ —; 2 



1 D f/2, 
2' 2"'"? 



(A4) 



(A5) 



For D = 2 one finds that the dependence of the energy 
on 7 is of the form E = ^"^ Eo + Q"^ f (U) , where Q is the 
total charge and /(J7) = (I/tt) d0 ln[l + {'^^-l)cos'^e]. 
Since we have to take a vanishing total charge, we are left 
with E = j'^Eq, so L = io7i which also conforms with 
eq.(A5) since F(l, i; 1; ^) = 7. Calculation for D = 1 
is also straightforward and the result is also given by 
eq.(A5): L = Lo7^- 



APPENDIX B: CALCULATION OF THE 
LAGRANGIAN FOR A POTENTIAL PARTICLE 

Here I calculate the Lagrangian for a spherical poten- 
tial particle. We need the middle term in eq.(27), which 
can be written as 



d''r'riif')q{j^) 



P_ 

27 



d^r'[Vfji/)Y, (Bl) 



)'U ■ Vijj['r{'r')] and fj that solves 



P_ 

27 

with qir") 

A'f] = g(r'). The integral is twice the electrostatic 
energy of the charge distribution q{r'), which is pro- 
duced by a polarized body with unidirectional polariza- 
tion -P(r') = —c~'^jtp['r{r')]U. We start from a spherical 
body having some V'(r) and divide it into thin spher- 
ical shells of radii and thicknesses dri of constant 
ip{ri). These are stretched into thin, concentric, nested 



homoeoids of minor axes and axes ratio 7 with con- 
stant polarization P{ri) along the major axis. We need 
the energy of this configuration. It is well known that the 
field inside an ellipsoid with uniform polarization along 
the major axis, is constant and proportional to the po- 
larization with the proportionality factor depending only 
on the axes ratio. In our case, for a constant if) wc write 
V'ry = d{j)P = —d{'-f)^c~'^ipU. (In the context of mag- 
netostatics ^(7) is called the demagnetizing factor.) This 
means that a thin homoeoid with uniform polarization P 
along its major axis produces a vanishing field inside it, 
and thus the interaction energy of two nested homoeoids 
such as ours vanishes. The Lagrangian produced by our 
stretched configuration is then the sum of the contribu- 
tions of the self energies of the individual thin homoeoids. 
Each thin homoeoid may be viewed as an infinitesimally 
thin dipole bilayer whose different elements thus do not 
interact with each other. Its energy is then the surface 
integral of the energy of a charged parallel-plate capac- 
itor E = (1/2)/ Saa^dS, where sigma is the surface 
density of the charge, and Sa is the thickness, both de- 
pendent on the polar angle 9' (relative to the major axis) 
at the position of the integration point q' on the ho- 
moeoid. Use instead as variable the polar angle 9 at 
point q on the spherical shell which was stretched into q'; 
we then have a = {U^ / c'^)'ii)cos9{cos'^9 + ^'^sin^9)^^/'^ , 
5a = jdri{cos'^9 + ^'^sin^Q)~^l'^ , dS = {cos^0 + 
j'^sm'^9y/'^nD-irf-^sin^-'^ede, where Qd is the D- 
dimensional solid angle Ojj = flo-i Jq sin^^^9d9. 

Thus the required contribution to L of the «th shell 
(for D >2) can be written as 



dLi = 



2 c' 



D-l 



cos'^9 + 72sm2| 



-de. 



(B2) 



Summing over the shells we replace ip^{ri)flD'ff' ^dr., 



The integral over 6 is 



r'^ cos^Osi 
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which can be read off formula 3.681.1 in GR to be 

r[(i^ - l)/2]r(3/2) D-l D + 2 _ 2. 
r[(L> + 2)/2] ^ ' 2 ' 2 ' ^ 

The factor in front can be shown to give ^Id/ D^d-i- I 
also use formula 9.131.1 in GR to transform 



Fa, 



D-l D + 2 



2, -2^n 3 D + 2 U\ 



2 ' 2 
and then 9.137.12 in GR to write 
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then adding the last term in eq.(27) I get finally 

L : 

with 



1 Z? [/2 



(B4) 



For = 1 it is straightforward to solve directly for 1] 
from cq.(24) and substitute in cq.(27). It turns out that 
eq.(B3) is still valid giving L = Lqj^. The same is true 
for D = 2 where we have L = Lq^. 

As a byproduct of the above calculation wc get the ex- 
pression for the D dimensional demagnetizing factor for 
a prolate ellipsoid magnetized along the symmetry axis. 
Consider a case where ip is constant inside the stretched 
ellipsoid in which V'ry = —d(j)c~'^^ipU . What 1 calcu- 
lated above is the quantity 

^ j d'^T f,{U ■ VV-) = / d^'r ^{U ■ Vt?) 
= d{j){j^p/2c^){U/cf J iP^d^r. (B5) 



This equals {p/2c'^){J iP'^d'^r)[F{l, f ; ^) - 1], as we 
found above. Comparing the two expressions one gets 

(B6) 
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3 this gives 
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rf(7) 



72-1 



7 



7 +(7^-1)^/^ 

2(72- 1)1/2"- \^-y-(-y2_ 1)1/2 



which reproduces the result found in [21]. 
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APPENDIX C: THE REST MASS OF A SOURCE 
PARTICLE 

The rest energy of a source particle is the energy differ- 
ence between two configurations, one of a uniform fluid 
at rest, the other likewise but with the source inserted. 
There are subtleties involved in the determination of this 
difference as the two configurations have infinite mass 
and energy. I adopt the following scheme: Consider a 
container of finite volume V much larger than that of 
the particle and filled with a static homogeneous fluid 
at the reference density p. Consider now a spherical 
source of vanishing total out-flux somewhere inside the 
container, in a steady state. Write the source density as 
s{r) = pA<J> = pV ■ u for some u = V<I>; and if the sup- 
port of s is within radius Rq of its center we deduce from 



Gauss theorem that $ and u vanish everywhere outside 
Rq. Writing the density as the continuity equation 

is V • [(p + C)^^ — P'w] = 0, with V and u radial from the 
center of the source; hence 



{l + C/p)-^u 



(CI) 



is exact and v vanishes everywhere outside the source. 
The Bernoulli equation is 



lv^ + h{p + 



constant = h{p). 



(C2) 



Since outside the source u = it follows that the density 
is constant there and equals p. 1 require this configura- 
tion with the source to have the same total fluid mass 
as the reference conflguration, so jy Qd^r = and this 
closes the set of algebraic equations that determines the 
configuration completely. The energy of the source con- 
figuration relative the reference one is then 



+ Cy+e{p + C)-e{p)]d''r. (C3) 



It is evident that neither the run of C and v inside the 
source nor the value of p, the constant density outside 
the source, depend on the position of the source inside 
the container, nor on the shape of the container (though 
its volume does enter). So far everything is exact. When 
the source may be considered a weak perturbation as in 
our case (C ^ p, w <C c) we can expand to lowest order 
in ^: Bernoulli's equation gives 



C(r)«C- 



pv"^ (f) 



(C4) 



where ( = p — p; the continuity equation gives tT ~ u 
V$; and the preserved-total-mass constraint gives 



p 



2c2y 



u^d 



(C5) 



where the integral is over the volume of the source. 
Clearly vanishes in the limit of infinite V. Writing 
e{p + Q — e{p) ~ h{p)(, the fact that ( integrates to zero 
means that there is no first order contribution to the in- 
trinsic energy difference. The second order contribution 
/ higher order in u^/c^. We are left with 

the dominant contribution 



J s 



- I ^sd^r 



fjsd r = Lq, 



(C6) 



which depends neither on the position of the source nor 
on the volume of the container or its shape, and we can 
identify it as the rest energy for large V. It is seen that 
Er = Eq/pSId, so we finally get 



D 



nio 
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One might also worry about the elastic energy of the 
container's wall, which is different in the reference 
configuration-whcrc the wall is subject to pressure p{p)- 
and the one with the source where the pressure on the 
wall is p{p) « p[p) + c^^. But this energy scales with the 
size of the container as oc V^^^"^ {A the area of the 
wall), so it can be neglected for large volumes. 



When more than one particle is present without over- 
lap they do not interact and the total energy is the sum 
of the rest masses. (Due to the non-linearity of the en- 
ergy there is interaction when particles overlap, but since 
$ = outside particles there is no long range interaction- 
similar to the case of two rigid spherical charge distribu- 
tions, each of vanishing total charge.) 



[1] W.G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 
[2] T. Jacobson, Phys. Rev. D 44, 1731 (1991). 

[3] M. Visser, gr-qc/9311028. 

[4] C. Barcelo, S. Liberati, and M. Visser, gr-qc/0505065. 
[5] W.G. Unmh, Phys. Rev. D 51, 2827 (1995). gr- 
qc/9409008. 

[6] M. Visser, Class. Quant. Grav. 15, 1767 (1998) gr- 
qc/9712010. 

[7] M. Milgrom, Astrophys. J. 270, 365 (1983). 
[8] M. Milgrom, Ann. Phys. 229, 384 (1994). astro- 
ph/9303012. 

[9] R.H. Sanders and S.S. McGaugh, Ann. Rev. Astron. As- 
trophys. 40, 263 (2002). astro-ph/0204521. 

[10] M. Visser and S. Weinfurtner, cond-mat/0409639. 

[11] M. Visser and S. Weinfurtner, Phys. Rev. D 72, 044020 
(2005). gr-qc/0506029. 

[12] G.E. Volovik The Universe in a Helium Droplet, Oxford 
University Press, England, (2003) 



[13] H. Lamb, Hydrodynamics, Dover, 6th edition (1932). 
[14] L.D. Landau and E.M. Lifshitz, Fluid Mechanics, Perga- 

mon Press, 2nd edition (1987). 
[15] M. Milgrom, J. Phys. A 35 (6), 1437 (2002). math- 

ph/0112040. 

[16] A.M.J. Schakel, Mod. Phys. Lett. B 10, 999 (1996). cond- 
mat/9607164. 

[17] M. Milgrom, Phys. Lett. A 190, 17 (1994). astro- 
ph/9401009. 

[18] R. Pciorls, More Surprises in Theoretical Physics, Prince- 
ton University Press (1991). 

[19] M. Stone, in "Analog models of General Relativity" 
Brazilian Genter for Research in Physics, Urea, Rio de 
Janeiro 16- 20th October 2000 cond-mat/0012316. 

[20] LS. Gradshteyn and I.M. Ryzhik, Tables of Integrals, Se- 
ries, and Products 5th edition, Academic Press (1994). 

[21] J. A. Osborn, Phys. Rev. 67, 351 (1945). 



